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ABSTRACT 



43 ■ 
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O ' Aims. We study the nucleosynthesis of several neutron capture elements (barium, europium, lanthanum, and yttrium) in local group dwarf 
spheroidal (dSph) galaxies and in the Milky Way by comparing the predictions of detailed chemical evolution models with the observed data. 



H 



i Methods. We compare the evolution of [Ba/Fe], [Eu/Fe], [La/Fe], [Y/Fe], [Ba/Y], [Ba/Eu], [Y/Eu], and [La/Eu] observed in dSph galaxies 
i and in our Galaxy with predictions of detailed chemical evolution models. The models for all dSph galaxies and for the Milky Way are able 
to reproduce several observational features of these galaxies, such as a series of abundance ratios and the stellar metallicities distributions. 
y—i • The Milky Way model adopts the two-infall scenario, whereas the most important features of the models for the dSph galaxies are the low 
^ ' star-formation rate and the occurrence of intense galactic winds. 

Results. We predict that the [s-r/Fe] ratios in dSphs are generally different than the corresponding ratios in the Milky Way, at the same [Fe/H] 
values. This is interpreted as a consequence of the time-delay model coupled with different star formation histories. In particular, the star- 
| formation is less efficient in dSphs than in our Galaxy and it is influenced by strong galactic winds. Our predictions are in very good agreement 
; with the available observational data. 
s ! ' Conclusions. The time-delay model for the galactic chemical enrichment coupled with different histories of star formation in different galaxies 
allow us to succesfully interpret the observed differences in the abundance ratios of s- and r- process elements, as well as of tt-elements in dSphs 
OO , and in the Milky Way. These differences strongly suggest that the main stellar populations of these galaxies could not have had a common origin 
' and, consequently, that the progenitors of local dSphs might not be the same objects as the building blocks of our Galaxy. 
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1 . Introduction that large galaxies formed through the central collapse of a gas 

cloud at high redshifts (Larson 1974; Sandage 1986; Arimoto 
Are the formation and evolution of the Milky Way (MW) linked & Yoshii 1987; Matteucci 1994; Calura, Matteucci & Menci 
to its neighbouring dwarf galaxies, in the sense that physical 2004), whereas the hierarchical cosmological scenario states 
processes (such as star-formation or galactic winds) occurring that large galaxies formed via the merger of minor blocks in 
in one type of galaxy influenced the physical processes occur- more recent epochs (Sommerville et al. 2001; Menci et al. 
ring in the other type? This a question that has been asked 2002). According to the second scenario, the dwarf spheroidal 
several times recently (Venn et al. 2004; Helmi et al. 2006; galaxies (dSph) of the local group could be the remains of the 
Geisler et al. 2007 and references therein) and a few clues blocks that gave rise to our galaxy. Therefore, our Galaxy and 
have emerged, but no definitive answers have given. The Local i oca l dSph galaxies should share some common evolutionary 
Group is an excellent laboratory for shedding some light on that features and properties, especially chemical signatures. A sim- 
matter, due to the possibility of studying the properties of dif- ii ar scenario was originally proposed by Searle & Zinn (1978) 
ferent types of galaxy in details. This includes large disk galax- f or the formation of the outer halo of the Milky Way. 
ies (such as the Milky Way) and dwarfs, in particular irregulars Qne possible approach t0 searching for evidence of this 
and spheroidals (Mateo 1998; Kunth & Ostlin 2000). joint evolution is t0 analyse the che mical history of stars in 
The knowledge of how local galaxies formed and evolved the Milky Way and in its satellite dSph galaxies (Venn et al. 
provides important clues to the formation of structures in the 2004). By adopting this approuch, Geisler et al. (2007) anal- 
universe and the whole subject of galaxy formation. The two ysed the scenario in which the Halo was formed by the same 
most popular scenarios for the assembly of galaxies contra- progenitors as the local dSph galaxies. They found no evidence 
diet each other, because the monolithical scenario proposes of this when comparing the chemical properties of stars in the 
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halo of the Milky Way and in local dSph galaxies, in particular 
the Sculptor galaxy. 

The evolution of chemical elements is normally studied 
by analysing the abundance ratios of certain key elements 
leading to a better understanding of galaxy evolution (Tinsley 
1980, Matteucci 1996). Comparing the observed patterns of 
these abundance ratios with predictions from chemical evolu- 
tion models helps constrain several parameters of the models, 
thereby providing clues to understanding important processes 
in the formation and evolution of galaxies and in the nucleosyn- 
thesis of chemical elements. 

Abundance ratios predicted by the models, unlike other ob- 
servables, strongly depend only on the adopted nucleosynthesis 
prescriptions, initial mass function (IMF), and star-formation 
rate (SFR). Consequently, comparing between models predic- 
tions and observations allows one to either constrain both the 
IMF and SFR of a galaxy if the nucleosynthesis of the elements 
analysed is well known or to obtain clues to the formation of 
the elements observed if the IMF and SFR are established. As 
a consequence, differences in the patterns of the same abun- 
dance ratios in different types of galaxies, on the other hand, 
reflect the particularities in the evolution of each type of system 
(Matteucci 1996). A comparison between abundance patterns 
in dSphs and the Milky Way therefore is necessary to ascertain 
wether these galaxies had a common origin or a completely 
different history of formation and evolution. 

In this work we analyse the trends of Ba and Eu, plus Y and 
La in the MW and in six local dSph galaxies (Carina, Draco, 
Sagittarius, Sculptor, Sextan, and Ursa Minor). By comparing 
the observed trends in these galaxies with chemical evolution 
models suited to adjusting their main features, we are able to 
impose constraints on the production of these elements and also 
to shed some light on the possible evolutionary connection be- 
tween the Milky Way and its neighbouring dSph galaxies. A 
comparison between the chemical evolution of the Milky Way 
and dSphs uses the standard model developed by Lanfranchi 
& Matteucci (2003, hereafter LM03). Previous papers (LM03, 
Lanfranchi & Matteucci 2004 (LM04)) already studied the 
[a/Fe] ratios in dSphs and predicted that they are generally 
lower than the corresponding ratios in the Milky Way at the 
same [Fe/H]. This is because the star formation in dSphs was 
assumed to have been much less efficient than in our Galaxy, 
thus producing the appearance of Fe from Type I a supernovae 
(SNe) at lower [Fe/H] values than in the Milky way. 

The paper is organised as follows. In Sect. 2 we present the 
observational data concerning the dSph galaxies and the Milky 
Way; in Sect. 3 we summarise the main characteristics of the 
adopted chemical evolution models for both types of galaxies, 
such as the star formation, the galactic winds and nucleosyn- 
thesis prescriptions; in Sect. 4 the predictions of our models 
are compared to observed abundance ratios and the results dis- 
cussed. Finally in Sect. 5 we draw some conclusions. All el- 
emental abundances are normalised to the solar values ([X/H] 
= log(X/H) - log(X/H) ) as measured by Grevesse & Sauval 
(1998). 



2. Data sample 

The abundance of several chemical species in stars of local 
dSph galaxies have become available in the past few years 
due to the advent of large telescopes and high-resolution spec- 
trographs (such as the Very Large Telescope - VLT - the 
FLAMES, and others). With this class of instruments several 
observational groups (Bonifacio et al. 2000, 2004; Shetrone, 
Cote & Sargent 2001; Shetrone et al. 2003; Venn et al. 2004; 
Sadakane et al. 2004; Fulbright, Rich & Castro, 2004; Geisler 
et al. 2005) have devoted an enormous effort to determining 
chemical abundances in red giant stars in a wide range of ages, 
making it possible to analyse the chemical enrichment history 
of these galaxies. From these observations we then gathered the 
data from the same galaxies in LM03 and LM04. Despite the 
relativly small number of data points, it is possible to compare 
the observed abundance ratios with the model predictions. We 
compared the evolution of [Ba/Fe], [Y/Fe], [La/Fe], [Eu/Fe], 
[Ba/Eu], [La/Eu], [Ba/Y], and [Y/Eu] observed in local dSph 
galaxies and in the Milky Way with the predictions of chemi- 
cal evolution models. 

Following Lanfranchi, Matteucci & Cescutti (2006a - 
LMC06a), we adopted the updated abundance values of Venn 
et al. (2004) for the dSph galaxies to have a more homoge- 
neous sample with data coming from different authors and to 
make a more proper comparison with the predictions of the 
models. In LMC06a, however, we only analysed the elements 
Ba and Eu. We extend the analysis to Yttrium and Lanthanum, 
for which there is a new set of data for Sagittarius. The abun- 
dances of La in Sagittarius, in contrast to Ba and Eu from 
Bonifacio et al. (2000), could have been underestimated due to 
the fact that HFS was not used by Sbordone et al. (2007 - pri- 
vate communication). In that case, the adopted values should 
be seen in the plots as upper limits. Like LMC06a, we ex- 
cluded some observed stars from the sample since they ex- 
hibit anomalous values for [Ba/H] or [Eu/H]. Two stars in Ursa 
Minor and two in Sculptor might have had their heavy ele- 
ment abundance affected by external factors. In Ursa Minor, 
the stars K and 199 (in Shetrone, Cote & Sargent 2001) exhibit 
heavy-element abundance ratios enhanced compared to those 
typical of other dSph stars: the K star has an abundance pat- 
tern dominated by the s-process and was classified as a carbon 
star, while Ursa Minor 199 is dominated by the r-process (see 
also Sadakane et al. 2004). In Sculptor, the stars Sculptor 982 
(Geisler et al. 2005) and Sculptor H-400 (Shetrone et al. 2003) 
exhibit enhanced heavy-element abundances. Both inhomoge- 
neous mixing of the SNe II (Shetrone et al. 2003) and enrich- 
ment by another star, which had already ejected its material 
in the ISM (Geisler et al. 2005), were invoked to explain the 
anomalous abundance patterns of the above stars. Since these 
stars do not exhibit any abundance pattern characterised only 
by the nucleosynthesis process occurring in the stellar interior, 
keeping them in the sample could lead to a meaningless com- 
parison with the model predictions and, as a consequence, to 
a misleading interpretation and to incorrect conclusions about 
the processes and the site of production of the heavy elements 
analysed. 
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For the Milky Way, we adopted the data taken from sev- 
eral works covering a wide range of metallicities from very 
metal-poor stars ([Fe/H] down to -4.0 dex) to stars with solar 
metallicities. The [Ba/Eu] ratio was collected from Francois et 
al. (2006), Burris et al. (2000), Fulbright (2000), Mashonkina 
& Gehren (2000, 2001), Koch & Edvardsson (2002), Honda et 
al. (2004), and Ishimaru et al. (2004). The other heavy element 
abundances were taken from Francois et al. (2006), Cowan et 
al. (2002), Burris et al. (2000), Johnson (2002), Pompeia et 
al. (2003), and Mc William & Rich (1994), in the case of lan- 
thanum, and from Francois et al. (2006), Burris et al. (2000), 
Mc William et al. (1995), Fulbright (2000, 2002), Mashonkina 
& Geheren (2001), Johnson (2002), Nissen & Schuster (1997), 
Prochaska et al. (2000), Gratton & Sneden (1994), Edvardsson 
et al. (1993), Stephens & Boesgaard (2002), and Honda et al. 
(2004) in the case of yttrium. 

3. Models 

In this work, we investigate the evolution of heavy elements 
in dSph galaxies and compare it to the evolution of the same 
elements in the Milky Way by means of detailed chemical evo- 
lution models, one suited to the dSph galaxies and the other 
to the Milky Way. The models for the dSph galaxies are the 
same as described in LM03 and LM04, whereas the model for 
the Milky Way is the one from Chiappini, Matteucci, Gratton 
(1997 - CMG97), as adopted by Cescutti et al. (2007). All mod- 
els adopt up-to-date nucleosynthesis yields for intermediate- 
mass stars (IMS) and supernovae of both types (type la and 
type II) as well as the effects of these objects in the energetics 
of the interstellar medium (ISM). The models allow the evolu- 
tion of the abundances of several chemical elements (such as 
H, He, C, N, O, Mg, Si, Ca, Fe, Ba, Eu, La, Y, Sr, Zr) to be 
followed in detail, starting from the matter reprocessed by stars 
and injected in the ISM through galactic winds and supernovae 
(SNe) explosions. 

The time evolution of the abundance in mass fraction of 
an element i in the gas of the galaxy is described by the basic 
equations given in Matteucci (1996) and Tinsley (1980). In par- 
ticular, the variation in time of the fractional mass of a chemical 
element, G„ is given by 

Gi = -<fiWi(t) + R,<t) + (Gdinf - (GiU, (1) 

where the gas mass in the form of an element ; normalised to 
a total fixed mass, M tot , is G,(f) = M g (t)Xi(t)/M tot and the to- 
tal fractional mass of gas present in the galaxy at the time t is 
G(f) = M g (t)/M tot . The abundance by mass of an element i is 
represented by the quantity X t (t) = Gi(t)/G(t), with the summa- 
tion over all elements in the gas mixture being equal to unity. 
The star formation rate (SFR), namely the fractional amount of 
gas turning into stars per unit time, is represented by if/(t), and 
specified by the star formation efficiency v, namely the inverse 
of the SF time-scale and expressed in Gyr 1 . Here, /?,■(?) repre- 
sents the returned fraction of matter in the form of an element i 
that the stars eject into the ISM through stellar winds and super- 
nova explosions, a term that contains all the prescriptions con- 
cerning the stellar yields and the supernova progenitor models. 
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The infall of external gas and galactic winds are accounted for 
by the terms (G,) 1M / and (G,) ouf , respectively. The type la SN 
progenitors are assumed to be white dwarfs in binary systems 
according to the formalism originally developed by Greggio & 
Renzini (1983) and Matteucci & Greggio (1986). 

3.1. The chemical evolution model fordSph galaxies 

In our scenario the dSph galaxies are formed through the in- 
fall of pristine gas until a mass of ~ 10 8 M© is accumulated. 
The star formation history (SFH) of each galaxy is given by 
the analysis of observed colour-magnitude diagrams (CMDs), 
which suggest a unique long (t ~ 3 to t ~ 8 Gyr) episode of 
SF in five systems, in particular, Draco, Sagittarius, Sculptor, 
Sextan, and Ursa Minor (Hernandez et al. 2000; Aparicio et 
al. 2001; Carrera et al. 2002; Dolphin et al. 2005). In the case 
of Carina, on the other hand, the SFH is characterised by four 
long (~ 2 Gyr of duration) episodes, as suggested in Rizzi et al. 
(2003), and successfully adopted in Lanfranchi, Matteucci & 
Cescutti (2006b - LMC06b). The episodes of SF, in our mod- 
els, are strongly affected by the occurrence of galactic winds. 
As soon as the wind starts, a large fraction of the gas reservoir 
that fuels the SF is removed from the galaxy, causing a sudden 
drop in the star formation rate. In that sense, the majority of 
galactic stars are formed before the galactic wind occurs, char- 
acterising this epoch as the most important period of galactic 
activity. 

The occurrence of the galactic wind depends on the binding 
energy of the galaxy and on the thermal energy of the gas, in 
the sense that when the second equates or it is larger than the 
first a galactic wind is deflagrated (see for example Matteucci 
& Tornambe 1987). The thermal energy strongly depends on 
the assumptions for the thermalisation efficiency of supernovae 
of both types and of stellar winds (Bradamante et al. 1998). In 
this work we assume, as standard values, the thermalisation ef- 
ficiencies suggested by Recchi et al. (2001), i.e. rj SNeIa = 1.0, 
rjSNeii = 0.03, 77s iv = 0.03 (for SNela, SNell and stellar winds, 
respectively). The binding energy, on the other hand, is strongly 
influenced by assumptions concerning the presence and distri- 
bution of dark matter (Matteucci 1992). A diffuse (RJR d =0.1, 
where R e is the effective radius of the galaxy and is the ra- 
dius of the dark matter core), but massive (Md ar k/M Lum = 10) 
dark halo has been assumed for each galaxy. This particular 
configuration allows the development of a galactic wind in 
these small systems without destroying them. 

The rate of gas loss G !VV , is proportional to the SFR, tf/{t), 
through a proportionality constant, w,, which is a free parame- 
ter describing the efficiency of the galactic wind: 

G iw = wnf/{t). (2) 

The values adopted for the efficiency of the wind (w ( ) are high 
(from 4 to 13) in order to reproduce the intense decrease ob- 
served in several abundance ratios, to remove a large fraction of 
the gas content of the galaxy, and to explain the observed stel- 
lar metallicity distributions (see Lanfranchi & Matteucci 2007 
for a more detailed discussion). The comparison with all these 
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Table 1. Models for dSph galaxies. 



galaxy 


v(Gyr~ l ) 


Wj n 


KCyr) 


d(Gyr) 


'aw (Gyr) 


IMF 


Carina 


0.15 


5 4 


0/2/7/9 


2/2/2/2 


0.53 


Salpeter 


Draco 


0.05 


4 1 





4 


1.97 


Sal peter 


Sagittarius 


3.0 


9 1 





13 


0.10 


Salpeter 


Sculptor 


0.2 


13 1 





7 


0.44 


Salpeter 


Sextan 


0.08 


9 1 





8 


0.78 


Salpeter 


Ursa Minor 


0.1 


10 1 





3 


0.43 


Salpeter 


Standard 


0.3 


10 1 





8 


0.33 


Salpeter 



constraints leaves little room for modifying the wind efficiency 
without changing the previous results much. 

The main assumptions of the dSph galaxy models are 

- the model is one zone with instantaneous and complete 
mixing of gas inside this zone; 

- no instantaneous recycling approximation is adopted, i.e. 
the stellar lifetimes are taken into account; 

- the evolution of several chemical elements (H, D, He, C, N, 
O, Mg, Si, S, Ca, Fe, Ba, Eu, La, Y, Sr and Zr) is followed 
in detail; 

- the nucleosynthesis prescriptions include the yields of 
Nomoto et al. (1997) for type la supernovae, Woosley & 
Weaver (1995) (with the corrections suggested by Francois 
et al. 2004) for massive stars (M > 10M G ), van den Hoek 
& Groenewegen (1997) for intermediate mass stars (IMS), 
and the ones described in Cescutti et al. (2006, 2007) and 
Busso et al. (2001) for Ba and Eu. 

The prescriptions for the SF (which follow a Schmidt law 
- Schmidt 1963), initial mass function (IMF - Salpeter 1955), 
infall, and galactic winds are the same as in LM03 and LM04. 
The main parameters adopted for the model of each galaxy, 
together with the predicted time for a galactic wind, tew, can be 
seen in Table 1, where v is the star-formation efficiency, w, the 
wind efficiency, n, t, and d are the number, time of occurrence, 
and duration of the SF episodes, respectively, t cw the time of 
the galactic wind. 

3.2. The chemical evolution model for the Milky Way 

We adopted the two-infall model (CMG97) for the evolution 
of the Galaxy. In this scenario, the Galaxy formed through two 
main episodes of infall: the first episode formed the halo and 
thick disk on a timescale of ~ 0.7 Gyr and the second one 
formed the thin disk on a timescale of ~ 7 Gyr for the so- 
lar vicinity. The formation of the thin disk is assumed to be a 
function of the galactocentric distance, leading to an inside-out 
build-up of the Galactic disk (see also Matteucci & Francois 
1989). Several concentrical rings, 2 kpc wide, with no ex- 
change of matter between them, simulate the Galactic disk. 
Two important characteristics of this model are (i) the almost 
independent evolution of the halo and the thin disk and (ii) the 
assumed threshold in the SF. When the gas density decreases 
below a critical value, in particular 7 M e pc~ 2 , the SF is halted 
(Kennicutt 1989, 1998; Martin & Kennicutt 2001). 

This model is able to reproduce a number of observational 
constraints of the solar vicinity and the whole Galaxy. The evo- 



lution of several elements as a function of [Fe/H] (such as C, N, 
a-elements, iron peak elements), several abundance ratios, the 
abundance gradients, and the G-dwarf metallicity distribution 
can be listed as the most important ones. In particular, Cescutti 
et al. (2007) used the same model to study the abundance gra- 
dients of several chemical elements including a-elements, iron 
peak elements, and heavy elements such as Ba, Eu, and La. 

The rate of mass accretion A(r,t) (the infall rate) is a func- 
tion of time and galactocentric distance and is given by 

A(r, t) = a(r)e-' /TH + b(r)e ( '- tm " )/TD(r) (3) 

where t max = \Gyr is the time for the maximum infall rate 
on the thin disk and th = 0.8Gyr and td are, respectively, 
the time scale for the formation of the halo thick-disk and the 
timescale of the thin disk. This last quantity is a function of the 
galactocentric distance: 

t d = 1 .033r(kpc) - 1 261 Gyr. (4) 

The coefficients a(r) and b(r) are constrained to reproduce 
the present-day total surface mass density as a function of 
galactocentric distance. In particular, b(r) is assumed to be dif- 
ferent from zero only for t > t max , where t max is the time of 
maximum infall on the thin disk (see Chiappini et al. 2003, for 
details). 

The adopted law for the SFR is substantially a Schmidt law 
with a dependence also on the total surface mass density. In 
particular, 

where the efficiency of the star formation process, v, is 
set to be lGyr 1 for the disk (f > IGyr). The total surface 
mass density is represented by 2(r, f), whereas E(r , f) is the 
total surface mass density at the solar position. The expres- 
sion G gas (r, t) represents the surface density normalised to the 
present time total surface mass density in the disk 2 D (r, teal), 
where teal = 13. 7 Gyr is the age assumed for the Milky Way 
and r Q = %kpc the solar galactocentric distance (Reid 1993). 
The exponent of the surface gas density, k, is set equal to 1.5. 
These choice of values for the parameters allow the model to 
fit the observational constraints very well, in particular in the 
solar vicinity. We recall that, below a critical threshold for the 
gas surface density (7M Q pc~ 2 for the thin disk and 4M Q pc~ 2 
for the halo phase), we assume that the star formation is halted. 
The adopted IMF is the one proposed by Scalo (1986). 
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Table 2. The stellar yields for La and Y in massive stars (r-process) in the case of primary origin. 
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12. I 9.00- 10- 8 I 4.00-10-' 

15. 3.00- 10-" 1.00- 10~ 8 

30. 1.00-10-' 1.00-1Q- 9 



4. Nucleosynthesis prescriptions 

The adopted nucleosynthesis prescriptions are the same as in 
previous works (see Lanfranchi, Matteucci & Cescutti 2006a 
- LMC06a) and in Cescutti et al. (2007). In particular, the pre- 
scriptions for the Fe and a-elements (namely O, Si, Ca, Mg) are 
those suggested in Francois et al. (2004). These authors anal- 
ysed several theoretical stellar yields by comparing of observed 
[el/Fe] vs. [Fe/H] trends with predictions for the Milky Way 
model in the solar neighbourhood and selected the sets of yields 
that fit the data best. For single low-intermediate mass stars, 
the stellar yields are those from van den Hoek & Groenewegen 
(1997) for the case of the mass loss parameter that varies with 
metallicity (see Chiappini et al. 2003, model5), whereas the 
yields of SNe la are taken from Iwamoto et al.(1997) with revi- 
sions in the case of Mg. The yields of SNe II that they suggest 
provides a best fit are the Woosley & Weaver (1995) ones. In 
this case, no modifications are required for the yields of Ca 
and Fe as computed for solar composition. For oxygen the best 
results are given by the Woosley & Weaver (1995) yields com- 
puted as functions of the metallicity. For the other elements, 
variations in the predicted yields are required to best fit the data 
(see Francois et al. 2004 for details). 

In the case of neutron process elements (s- and r-process), 
we adopted the yields of Cescutti et al. (2006) and LMC06a 
for Ba and Eu and the yields of Cescutti et al. (2007) for La. 
In these works, the yields of s-process elements (in particular 
La, Ba, and Y) are those predicted by Busso et al. (2001) in 
the mass range 1.5 - 3M Q , with an extension to the mass range 
1.5 - 1M Q . The r-production of Ba and Eu were obtained by 
following the prescriptions of model 1 of Cescutti et al. (2006). 
These yields are empirical and were chosen to provide the best 
fit to the abundances of Ba and Eu in low-metallicity stars in the 
Milky Way, as measured by Francois et al. (2007). The same 
set of yields were also adopted by LMC06a in the compari- 
son between the predictions of the dSph galaxy models with 
the abundances of Ba and Eu in red giant stars in local dSph 
resulting in a good agreement. 

Barium is assumed to be mainly produced by s-process 
in low and intermediate-mass stars (LIMS) in the mass range 
1 - 3M Q with a minor fraction also produced as an r-process 
element in massive stars (12 - 30 M ). Europium on the other 
hand is considered to be a pure r-process element produced in 
massive stars in the same range of masses. For La, we followed 
the prescriptions of Cescutti et al. (2007), who adopted the 
same prescriptions of Ba: a major contribution from s-process 
in LIMS and a minor production as r-process in massive stars 
in the same mass range as Ba. The yields for the s-process are 
taken from Busso et al. (2001), whereas those from r-process 
are the ones from Cescutti et al. (2007) (see Table 2). 



5. Results 

It has been suggested that Local Group dSph galaxies could 
be the remaining blocks of the small systems that gave rise 
to our Galaxy (White & Rees 1978; Navarro, Frenk & White 
1997). From the chemical evolution point of view, one can 
analyse such possible connection by comparing the abundance 
and abundance ratios of several elements and the metal-poor 
tail of the stellar metallicity distributions observed in the Milky 
Way and in local dSphs (Venn et al. 2004; Ripamonti et al. 
2007; Helmi et al. 2006, Geisler et al. 2007). It is expected, 
for example, that the s- and r-process ratios in dSph galaxies 
might be different from the same ratios in stars with compa- 
rable metallicities in our Galaxy, due to different SFHs. In the 
dSphs the SFR is much lower, so the contribution from the s- 
process should appear at lower metallicities when compared 
to the Milky Way, because of the slower increase in [Fe/H] 
with time. Venn et al. (2004) examined the abundance and 
abundance ratios of several s-process and r-process elements in 
dSph and in the Milky Way stars and claim that no significant 
stellar component of our Galaxy appears to have been formed 
from the dSphs. Helmi et al. (2006) uses a completely differ- 
ent approach to reach the same conclusion. They compared the 
metal-poor tail of the stellar metallicity distributions of four 
dSph galaxies to the one of the Galactic halo and found signifi- 
cant differences. The differences, according to the authors, sug- 
gest that the progenitors of local dSph galaxies might have been 
fundamentally different from the building blocks that formed 
our Galaxy. 

We intend not only to examine such a subject but also to 
investigate the different evolutions of heavy elements in these 
galaxies in order to constraint the formation and evolution of 
these elements and the evolution of these dwarf galaxies. By 
means of the detailed chemical evolution models described 
above, we are able to predict the evolution of La, Ba, Y, and 
Eu. The predictions for the dSphs allow us to fill the gaps in 
the observations, given the reduced number of available data. 
At high or very low metallicities, for example, there is almost 
no observed data for the dSphs, which prevents a comparison 
to extend to the whole range of metallicities observed in our 
Galaxy. 

For the Milky Way we applied the two-infall model 
(CMG97) as adopted by Cescutti et al. (2006). In the case of 
the dSph, galaxies we adopted the standard model of LM03. 
This model was constructed to reproduce the observed [a/Fe] 
ratios of a sample of dSph galaxies (Carina, Draco, Fornax, 
Leol, Leo II, Sagittarius, Sculptor, Sextan, and Ursa Minor) as 
a whole, without separating each individual galaxy. We con- 
sider that the general trends of the heavy elements in the dSph 
galaxies are the same for all individual galaxies, in spite of the 
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Fig. 1. The evolution of [La/Fe] and [Y/Fe] predicted by the 
dSph model for all six galaxies. 

observed differences caused by the particularities in the evolu- 
tion of each galaxy (in particular, different SF efficiencies and 
galactic wind efficiencies). By doing that, we intend to high- 
light the differences between the pattern of these elements in 
the Milky Way and its spheroidal satellites. There are, conse- 
quently, a few points that lie above or below the model's predic- 
tion, depending on the characteristics of the galaxy (SF history, 
SFR, galactic wind efficiency). Since the aim of the standard 
model is to draw a general trend for all analysed dSph galaxies 
, small offsets between the predictions from the standard model 
and the data of individual galaxies should not be seriously taken 
into account. The goal is to highlight the differences between 
the Milky Way and the dSph galaxies as a whole. 

We should mention, however, that Sagittarius dSph galaxy 
should be seen as an exception, not only for the particular 
chemical characteristics (elemental abundances, abundance ra- 
tios - Sbordone 2006 - star formation efficiency, and metallic - 
ity distribution - LM04), but also for the fact that this galaxy 
is clearly interacting with our Galaxy (Mateo et al. 1996; 
Majewski et al. 2003.) The other dSph galaxies exhibit simi- 
lar chemical properties (similar values for the chemical abun- 
dances, abundance ratios, SF efficiencies - LM04) with small 
differences due to their particular formation and evolution. The 
same properties in Sagittarius, on the other hand, are more like 
our Galaxy than like other dwarfs (see also Geisler et al. 2007). 
Consequently, when comparing local dSph galaxies with the 
Milky Way galaxy, we decided to exclude Sagittarius. 

In Figure 1 we show the predictions of La and Y in the 6 
dSphs. There is no available data for La in Draco, Ursa Minor, 
and Sextan, so we only show the model predictions in these 
cases. The predicted and observed trends of these two elements 
can be explained by the s- and r-process contributions from 
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Fig. 2. The evolution of [La/Fe] and [La/Eu] predicted by the 
Sagittarius model with La yields in LIMS increased by a factor 
of 2 and 4 (dashed and solid lines, respectively). We ruled out 
the model with the La yields increased by a factor of 4, even 
though it provides a better fit to the dSph observations, because 
it cannot reproduce the Milky Way data. 



stars in different mass ranges and by the occurrence of strong 
galactic winds. At low metallicities, the r-production in mas- 
sive stars plays a dominant role by being responsible for the 
plateau at almost solar values for [La/Fe] and [Y/Fe]. At inter- 
mediate to high metallicities ([Fe/H] > -1.8 dex), the injection 
of these elements in the ISM by LIMS stars becomes impor- 
tant in the dSph galaxies. Besides that, at [Fe/H] — 1.6 dex, a 
galactic wind is developed, thereby removing a large fraction 
of the gas content of the galaxy and reducing the SFR consid- 
erably. With a low SFR, the contribution from massive stars is 
almost negligible, whereas the one from LIMS remains impor- 
tant. These combined effects result in a decrease in the abun- 
dance ratios. 

The difference between the trends among these elements 
are related to details of their nucleosynthesis, in particular the 
different amounts of the fractional mass produced by each pro- 
cess for each element. The fraction of Y produced in LIMS by 
s-process is lower than the one of La. Consequently, one can 
observe a very similar trend of [La/Fe] and [Y/Fe] in all galax- 
ies, but with a more intense decrease in [Y/Fe] after the occur- 
rence of galactic winds, due to its lower production in LIMS. 

In all six dSph galaxies we are able to satisfactorily repro- 
duce the observed data for all elements, maybe with the ex- 
ception of La in Sagittarius. In this case, there is a discrepancy 
between model predictions and observations, in the sense that 
we under-predict the values of [La/Fe] at high metallicities (see 
Figure 1). It seems, in fact, that there is a trend in the observed 
data of increasing [La/Fe] with increasing metallicity, whereas 
in the model predictions one notices a slight decrease. The de- 
crease is a result of the effects of galactic winds on the SFR, 
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Fig. 3. The evolution of Ba, La, Y, and Eu in the Milky Way and local dSph galaxies. The thick line corresponds to the predictions 
of the Milky Way model and the thin line the predictions of the standard model for the dSph galaxies. 



as mentioned before. Could this discrepancy be solved with a 
change in the value of the galactic wind efficiency? Probably 
not. If we change the efficiency of the wind for Sagittarius, 
we would lose the agreemeent with other abundance ratios, not 
only with the [alpha/Fe] ratios but also with other heavy ele- 
ments. From Figure 1, it is clear that if we increase the wind 
efficiency (it would be necessary to match La data), the agree- 
ment with [Y/Fe] would be lost. In that case, we should look 
for another explanation for this discrepancy. In that metallicity 
range, the production of La is almost totally dominated by the 
s-process in LIMS. Consequently, an under-prediction of La in 
LIMS could be the source of the discrepancy. We adopted the 
yields of Busso et al. (2001) for La in this mass range, which 
could be increased by hand to get a better fit to the observa- 
tions. In that case, however, we would face a serious problem. 
The same set of yields when adopted in the model of the Milky 
Way allows very good agreement with the observed data, in- 
cluding the solar value (see Cescutti et al. 2006). Consequently, 
if we change the yields considerably (for a factor of 2 or more) 
we would lose the agreement with the Milky Way data and the 
solar value. 

If the yields cannot be considerably modified, the other 
source of discrepancy could reside in the data. The data for La 
in Sagittarius at high metallicity comes entirely from the work 
of Sbordone et al. (2006). In such work the authors did not use 
hyper-fine splitting (HFS) in the calculation of La abundances, 
which could cause an over-prediction of the suggested values, 
but not as high as the difference between the observed data and 
the model predictions (Sbordone, private communication). In 
that sense, the combined effects of overestimating the observed 



values with the possible under-prediction of La yields in LIMS 
could be the cause of the discrepancy between the observed 
data and model predictions. In order to examine this scenario, 
we ran two alternative models in which the La yields in LIMS 
were increased by a factor of two and four (Figure 2). One can 
notice that an increase in the yields of La in LIMS by a factor of 
two is not enough to reach the observed high values, which are 
only fitted by the model with yields four times higher than the 
ones of Busso et al. (2001). Since the agreement obtained with 
the solar neighbourhood and the solar values is lost if we adopt 
yields increased more than two times, then it is highly probable 
that there is also a overestimation of the observed values. 

The comparison between the predictions of the models and 
the observed values for neutron capture elements strongly sug- 
gests that La and Y trends can be explained in all dSph galax- 
ies (with maybe the exception of Sagittarius - see discussion 
above) by models that adopt a scenario in which both these ele- 
ments are produced by the s-process in LIMS, in the mass range 
M = 1 - 3M Q , with a minor fraction produced by the r-process 
in massive stars (M = 12 - 30M G ), following closely the pro- 
duction of Ba. Strong galactic winds also play a fundamental 
role in explaining the behaviour of La and Y at intermediate to 
high metallicities, as in the case of Ba, Eu, and a elements. 

In Figure 3, we show the predictions of the models of both 
dSphs and the Milky Way for Ba, La, Y, and Eu along with 
the observed data for the studied systems (small open circles: 
Milky Way stars, large symbols: dSph galaxies). This model 
was constructed to reproduce the observed [a/Fe] ratios of a 
sample of dSph galaxies (Carina, Draco, Fornax, Leol, Leo 
II, Sagittarius, Sculptor, Sextan, and Ursa Minor) as a whole, 
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Fig. 4. The evolution of [Ba/Y], [Ba/Eu], [La/Eu], and [Y/Eu] in the Milky Way and local dSph galaxies. The thick line corre- 
sponds to the predictions of the Milky Way model and the thin line to the predictions of the dSph galaxies standard model. 



without separating each individual galaxy, so this model repre- 
sents a typical dSph galaxy. 

At a first glance, it seems that the data from the two types of 
galaxies overlap at intermediate metallicities, which could sug- 
gest a common origin. A more careful analysis, however, re- 
veals that there are marked differences. In both types of galax- 
ies, the observed values show a large dispersion at low metallic- 
ities, which makes the interpretation more complex. This wide 
spread is normally attributed to an inhomogeneous mixing; i.e., 
the abundance pattern of r-process elements of a newly-formed 
star was contaminated by SNe II explosion in the vicinity. Even 
with the dispersion, one notices that, while the abundances of 
La, Ba, and Eu are similar in the Milky Way and in the dSph 
galaxies, [Y/Fe] exhibits significant lower values in the major- 
ity of the stars in dSph galaxies when compared to the ones of 
the Milky Way. This offset in [Y/Fe] suggests that the enrich- 
ment from r- and s-process elements in the two types of galax- 
ies might have been different, as already suggested by Venn et 
al. (2004). Besides that, the predictions of the models suggest 
meaningful differences at low and, especially, at intermediate- 
to-high metallicities. 

At low metallicities ([Fe/H] — 2.5 dex), there are several 
stars with low abundances of [La/Fe], [Ba/Fe] and [Eu/Fe] in 
the Milky Way, whereas these elements exhibit solar values in 
the dSph. There are a few stars in the dSph that show low val- 
ues of [Ba/Fe], similar to the ones in the Milky Way, but they 
could not be reproduced by the model with the same nucleosyn- 
thesis prescriptions adopted in the model of our Galaxy. Since 
these stars belong to different galaxies (Ursa Minor, Draco, and 
Sextans) and come from different works (Venn et al. 2004, who 



homogenized the atomic data for spectral lines, and Sadakane 
et al. 2004), it is not likely that the low [Ba/Fe] values are 
representative of one specific galaxy or that they are caused 
by a systematic error. In that way, these low values, if con- 
firmed by further observations, could only be reproduced if a 
different production of Ba is adopted. In such a case, the r- 
production of Ba in massive stars should be ignored and the s- 
production should be extended to stars with masses up to 4M© 
(see LMC06a). In either case, the metal-poor stars that enriched 
the ISM of the dwarfs and the MW could not be the same, sug- 
gesting that these metal-poor stellar populations are different 
(see also Helmi et al. 2006). If one considers the predictions 
of the models, the two populations are characterised by differ- 
ent abundance patterns, but if one considers the stars observed, 
then the production of heavy elements is not the same. It should 
be noticed, however, that the numbers of metal-poor stars ob- 
served in local dSphs are still small and further observations 
are required to confirm the low values of [Ba/Fe]. These two 
facts suggest that if metal-poor stars in local dSph galaxies do 
exist with low abundances of neutron capture elements, similar 
to the values observed in our Galaxy, then a different mecha- 
nism for the production of r- and s-process elements in dSph 
galaxies should be invoked. 

From intermediate-to-high metallicities, there are almost 
no stars observed in local dSphs, but the predicted trends are 
quite different from the one in our Galaxy. In the dSphs, the 
abundances of neutron capture elements (as all other elements) 
are strongly influenced by the development of strong galactic 
winds. This decreases the abundance of r-products but not as 
much as the s-products, giving rise to very low [Eu/Fe]. In the 
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Fig. 5. Comparison between the evolution of [Ba/Y], [Ba/Eu], [La/Eu], and [Y/Eu] predicted by the Milky Way model and the 
Draco and Sculptor models. The thick solid line corresponds to the predictions of the Milky Way model, the thin solid line to the 
Sculptor model, and the thin dashed line to the Draco model. The data from Sculptor are represented by solid squares, the data 
from Draco by solid triangles, and all the other dSphs by the solid circles. 



Milky Way, the feature that plays a major role in the abundance 
patterns of neutron capture elements at high metallicity is the 
s-production of these elements by LIMS. Since these stars have 
a long lifetime (up to a few Gyrs), they only contribute to the 
enrichment of the ISM when [Fe/H] has had time to increase 
to intermediate or high values, depending on the SFR. In the 
Milky Way, the injection of s-products from LIMS marked by 
the change in the slope of the [Ba/Eu] ratio both in the data 
and in models predictions) starts becoming important at around 
[Fe/H] — 1.0 dex, whereas in local dSphs, it happens at lower 
metallicities, around [Fe/H] ~ -1.7 dex. The reason for this dif- 
ference in the model can be explained by the fact that, even 
though Ba is injected into the ISM on the same timescale in 
both types of galaxies, the dSphs are less enriched in [Fe/H] 
than the ISM in the Milky Way at a particular time, owing to 
the less efficient star formation assumed in the dSphs compared 
to to the Galaxy. This is an important result that derives directly 
from the time-delay model (delayed Fe production relative to 
alpha-elements) coupled with different SFHs (see Matteucci 
2003). The changes in the abundance patterns of heavy ele- 
ments at different metallicities in the Milky Way and in the 
dSph, due to the injection of fresh elements by LIMS into the 
ISM, also strongly suggest that the two metal-poor populations 
on the different types of galaxies are not the same and do not 
have a common origin. 

In Figure 4, we show the predictions of the models of both 
galaxies for [Ba/Y], [Ba/Eu], [La/Eu], and [Y/Eu], along with 
the observed data for the studied systems. The analysis of abun- 



dance ratios between heavy elements strengthens the above 
scenario. At high metallicities one can notice an increasing 
trend of [Ba/Eu], [La/Eu], and [Y/Eu] in dSphs much more in- 
tense than the ones observed in the Milky Way. In our Galaxy, 
the increase in [Ba/Eu] is normally attributed to contributions 
of the s-process from LISM, whereas in dSph the interruption 
of the enrichment of r-process elements from massive stars due 
to the galactic winds plays a major role. The ratio [Ba/Y] re- 
flects the production of light and heavy s-process elements, 
since Ba is almost totally produced by s-process and the frac- 
tion of Y produced by s-process is smaller. In the Milky Way, 
the [Ba/Y] values are always lower than in the dSph galaxies, 
indicating different productions of these two elements in both 
types of galaxies. In the same way, several stars in the dSph 
galaxies exhibit lower [Y/Eu] abundance ratios when compared 
to stars with the same metallicity in our Galaxy. The nucleosyn- 
thesis of r- and s-process elements also seems to have followed 
different tracks in each type of galaxy, as is seen in the analy- 
sis of [Y/Eu]. From low-to-intermediate metallicities, [Y/Eu] is 
lower in the dSph galaxies than in the MW, suggesting different 
amounts of r- and s-process elements in each type of galaxy. 

In Figure 4, one sees that the prediction of the standard 
model for dSph for [Ba/Y] does not reach the highest values 
observed. As mentioned before, this model was adjusted to fit 
the data of all dSph galaxies at the same time, leaving behind 
the differences (in particular v,- and wf) among the systems. The 
stars with the highest values of [Ba/Y] belong to Draco, galaxy 
that is characterised by a much lower SF efficiency (v* = 0.05) 
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and by a lower wind efficiency (w, = 4). These values allow 
the specific model for Draco to reproduce the highest [Ba/Y], 
which are not fitted by the standard model (Figure 5). In Figure 

5, we also show the predictions of the Sculptor dSph galaxy. By 
comparing the models for individual galaxies with the model 
for the Milky Way, it can be seen the differences are even more 
pronounced. In particular, the differences between the predic- 
tions from both Draco and Sculptor models and the Milky Way 
model are even larger (when compared to the standard model) 
for [Ba/Y], [Ba/Eu], and [La/Eu] at high metallicities. 

Several differences in the evolutionary tracks of Ba, La, 
Eu, and Y, in the MW and local dSph galaxies are noticed in 
the comparison of their evolution from both the observational 
point of view and model predictions. The much lower SFR 
and the intense galactic winds in local dSph galaxies are the 
main features responsibles for the marked differences at low 
([Fe/H] < -2.0 dex) and high ([Fe/H > -1.6 dex) metallicities. 
Consequently, we suggest that the very old, metal-poor, stellar 
populations present in these galaxies are not the same and do 
not share a common origin. 

6. Summary 

We analysed the evolution of several r- and s-process elements 
(in particular barium, europium, lanthanum, and yttrium) in 
Local Group dwarf spheroidal (dSph) galaxies and in the Milky 
Way by comparing observed data in a wide range of metallic- 
ities with predictions of detailed chemical evolution models. 
The adopted models for all dSph galaxies and for the Milky 
Way are the ones described in LM03 and Cescutti et al. (2006), 
respectively. The Milky Way model adopts the two-infall sce- 
nario, whereas the most important features of the models for 
the dSph galaxies are the low star formation rate and intense 
galactic winds. By assuming that Eu is a pure r-process ele- 
ments synthesised in massive stars (M — 12 - 30M o ), whereas 
Ba, La, and Y are produced mainly by s-process in low and 
intermediate-mass stars (LIMS) in the mass range M = 1 -3M 
with a minor contribution from r-process in massive stars, we 
are able to reproduce the observed abundance trends of all 
abundance ratios very well. At low metallicities the trends are 
defined by the contribution from the r-process in massive stars 
in both galaxies, whereas at high metallicities the s-process in 
LIMS are the main contributor to the enrichment of these el- 
ements in our Galaxy and in the dSph galaxies. In these last 
ones, however, strong galactic winds play a fundamental role 
at [Fe/H] > -1.6 dex. When comparing the data and predictions 
of the Milky Way model with the standard model for the dSph 
galaxies, significant differences both in the observed trend and 
in models predictions at all metallicities are seen. These differ- 
ences strongly suggest that the metal-poor stellar populations 
are not the same in dSphs and the Milky Way. 
The main conclusions are the following: 

- The observed trends of [La/Fe] and [Y/Fe] are reproduced 
very well by the models of all six dSph galaxies and by the 
Milky Way model, suggesting that our assumed nucleosyn- 
thesis prescriptions are reasonable. In particular, we assume 
that La and Y are are mainly s-process elements produced 



in LIMS (M = 1 - 3M ) with a minor fraction produced by 
r-process in massive stars in the range M — 12 - 3OM . 

- The Sagittarius model seems to underpredict the values of 
[La/Fe] at high metallicities when compared to observa- 
tions. One possible source of this disagreement could be 
an underestimation of La yields in LIMS by Busso et al. 
(2001). A model with an increased yield of La by a factor 
of two can reproduce the data very well if they are overesti- 
mated by no more than ~ 0.2 dex. An increase greater than 
a factor of two would, however, cause a disagreement with 
the Milky Way data, so it is not likely that this is the only 
source of disagreement. 

- The predictions for [Ba/Fe], [Eu/Fe], [La/Fe], and [Y/Fe] 
by the dSph standard model developed by LM03 and the 
Milky Way model reveals several differences, as suggested 
by observed data. At low metallicities ([Fe/H] < -2.0 dex), 
the models of dSph galaxies predict higher values for these 
ratios due to the lower SFR, whereas) they are always lower 
at high metallicities ([Fe/H] > -1.4 dex, due to the effects 
of galactic winds on the SFR. At intermediate metallicities, 
[Ba/Fe], [Eu/Fe], and [La/Fe] seem to be similar in both 
types of galaxies but [Y/Fe] is lower in the dSph galaxies. 

- The heavy element abundance ratios ([Y/Eu], [Ba/Eu], 
[La/Eu], and [Ba/Y]) also exhibit different trends in our 
Galaxy and its neighbour dwarfs: at high metallicities, the 
values in dSph galaxies are always higher, whereas at inter- 
mediate metallicities [Y/Eu] is lower and [Ba/Y] is higher 
in the dSphs. 

- There are a few dSph stars with very low [Ba/Fe] at low 
metallicities, similar to what is observed in stars with com- 
parable metallicities in our Galaxy. If these values are in- 
deed real, a different nucleosynthesis of Ba in these galax- 
ies is required. In such a case, the r-production of Ba in 
massive stars should be suppressed and the s-production 
should be extended to stars with M = 4M G (see LMC06a). 

- The meaningful differences between the trends of heavy 
elements in local dSph galaxies and our Galaxy strongly 
suggest that the stellar populations of these two types of 
galaxies are different and might not have shared a common 
evolution. In fact, the evolution in the dSph galaxies must 
have proceeded at a much lower rate than in the Milky Way 
and been strongly influenced by intense galactic winds. In 
that sense, it is unlikly that the progenitors of local dSph 
galaxies and the building blocks of the Milky Way are the 
same objects. 
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